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ABSTRACT 

We obtained six observations of PSR J1741-2054 using the Chandra ACIS-S detector totaling ~300 ks. 

By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray 
point sources in the field of view, we have measured the pulsar proper motion at p = 109 ± lOmas yr 1 in a 
direction consistent with the symmetry axis of the observed Ha nebula. We investigated the inferred past 
trajectory of the pulsar but find no compelling association with OB associations in which the progenitor may 
have originated. We confirm previous measurements of the pulsar spectrum as an absorbed power law with 
photon index T=2.68±0.04, plus a blackbody with an emission radius of (4.5 4 2 5M0.38 km, for a DM-estimated 
distance of ().38<7o.3x kpc and a temperature of 61.7 ±3.0 eV. Emission from the compact nebula is well 
described by an absorbed power law model with a photon index of T = 1.67 ±0.06, while the diffuse emission 
seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also applied 
image deconvolution techniques to search for small-scale structures in the immediate vicinity of the pulsar, but 
found no conclusive evidence for such structures. 

Subject headings: pulsars: individual (PSR J1741-2054) - X-rays: general 


1. INTRODUCTION 

PSR J1741-2054 (J1741) is one of the closest middle-aged 
(t c = 390 kyr) pulsars known. It has a period of P = 413 ms 
and was first discovered in 7-rays using the Large Area Tele¬ 
scope (LAT) on the Fermi Gamma-ray Space Telescope by a 
blind search f or peri odic 7-ray pulsations from Fermi- LAT 
point sources (lAhdo!l2009b . It was subsequently detected in 
archival Parkes_radio data and observed using the Green Bank 
Telescope (ICamilo et al.ll2009l) . The pulsar has a spin down 
energy loss rate of E = 9.5 x 10 33 erg s _1 which is moderately 
low compared to those of other 7-ray pulsars. The pulsar has a 
very small dispersion measure (DM) = 4.7 pc cm' 1 2 3 and a mag¬ 
netic field of 2.7 x 10 12 G. Using the NE2001 Galactic electron 
density model, the low DM implies a distance of 0.38 kpc 
(tCordes & Lazioll2002i) . At this distance, its measured radio 
flux at 1400 MHz (S ~ 160 pjy) makes it one of the least lumi¬ 
nous radio pulsars known. Its 7-ray pulsations lag behind its 
radio pulsations by 6=0.29P, implying that our line-of-sight 
tangentially cuts the 7-ray cone, while nearly missing the ra¬ 
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dio beam (iCamilo et all 120091) . This makes J1741 a transi¬ 
tional object between a classical radio/7-ray loud pulsar such 
as Vela_and the radio- quie t Geminga-type pulsars. Interest¬ 
ing I y, iRomani et al.l(1201 Ol ) detected a 20" long, non-radiative 
Ha- bow shock nebula around the pulsar. Modeling of the bow 
shock suggested that the pulsar is traveling with a velocity of 
~ 150 km s“\ while the observation of negative radial veloc¬ 
ities from both sides of the nebula imply that the velocity is 
directed out of the plane of the sky at an angle of 15° ± 10°. 

Using a short Chandra ACJX-S observation (observation 
ID (ObsID): 11251), IRomani et all (1201 Oh detected an X- 
ray pulsar wind nebula (PWN) within this Ha nebula and a 
long (> 2 arcmin ) X-ray tr ail at an angle of 45° ±5° East 
from North. IRomani et al.l (|2010[) also suggested that there 
are asymmetries in the small scale structure surrounding the 
pulsar, which they associate with a com pact 2.5" equatorial 
toroid al structure. Ma relli et al.l (120141) and iKarnova et al.l 
(120141) performed a spectral analysis of the pulsar emission 
using XMM-Newton and Chandra and determined that a two- 
component (blackbody plus power-law) model is required to 
obtain satisfactory spectral fits. 

In this paper, we use ~300 ks of Chandra data of J1741 
that were obtained as part of the Cycle 14 Chandra Vision¬ 
ary Project “A Legacy Study of the Relativistic Shocks of 
PWNe”, plus a ~ 49ks archival observation, to constrain the 
pulsar motion, and the geometry of the PWN outflow. We 
discuss our approach to image registration and proper motion 
measurement in Section 3, followed by a discussion in Sec¬ 
tion 4 of our image deconvolution efforts to search for small- 
scale structure around the pulsar. In Section 5 we discuss the 
results of our spectral fits for the pulsar emission and that of 
the extended PWN, and compare these with previous results. 
In Section 6 we discuss the implications of the proper motion 
measurements, including comparisons with the observed Ha 
nebula surrounding J1741, and discuss the lack of evidence of 
for synchrotron cooling in the PWN trail. Our conclusions are 
summarized in Section 7. 






























2 



06 04 02 17 : 42:00 58 56 54 41:52 


Right ascension 

Figure 1. Merged Chandra exposure-corrected 0.3-5.0 keV ACIS-S image of the extended emission around J1741-2051. This was produced using repro ject_obs 
and flux_obs and incorporates all available Chandra observations. The image is smoothed using a Gaussian of width 3" and plotted on the logarithmic scale. 
The reference sources numbered 1 through 13 were used for relative astrometry, while the pulsar is labeled as P. Other sources seen in the image are not suitable 
for astrometry because, being variable, they were not significantly detected in both the archival observation and one of the new epoch observations. We included 
a cutout region around the pulsar to show the point source. The magenta lines show our derived proper motion (solid) and its uncertainty (dashed), traced back 
from the pulsar position. The background region used for spectral analysis is shown as the magenta dashed circle. 


2. X-RAY DATA ANALYSIS 

We obtained 282 ks of new Chandra ACIS-S exposure time 
of J1741. TableQ]lists the parameters for each of the six new 
observations that we obtained, including the archival observa¬ 
tion. To reduce pileup, the CCDs were operated in half-frame 
VFAINT mode so that events were read out every 1.7s. The 
pulsar was placed near the optimum focus on the backside il¬ 
luminated S3 chip. In addition to our six observations, there 
is a 48.8 ks archival observation (ObsID: 11251), which was 
taken on 2010 May 21. Each new observation has a roll angle 
similar to the roll angle of the 2010 epoch (~ 90°), except for 
ObsID: 15544. This observation has a roll angle of ~ 260° as 
the nominal roll angle of Chandra is rotated by 180° during 
the time of the year this observation was completed. The data 
were analysed with CIAO 4.6.2 after all observations were re¬ 
processed using the CALDB 4.5.9. No flaring occurred in any 
of the observations so the full exposure times were used. 

Using all seven observations we produced a merged, 
exposure-corrected image of J1741 by reprojecting the new 
observations to a common tangent plane based on the WCS 


information of ObsID: 11251 (CIAO task: reproject_obs) 
and combined all reprojected observations into an exposure 
corrected image using the CIAO task flux_obs. The merged 
ACIS-S image of the extended emission around the pulsar, 
smoothed with a 3" Gaussian", is shown in Figure |T| The 
pulsar point source lies at the apex of the diffuse X-ray emis¬ 
sion, while a diffuse, faint X-ray trail extending ~ 1.9' is seen 
towards the northeast of the pulsar. 

3. REGISTRATION AND THE PROPER MOTION 

To constrain the proper motion of the pulsar, we registered 
each of the new Chandra images to the archival image using 
field point sources that were identified using the CIAO tool 
wavdetect. We selected sources with a detection significance 
of > 3cr that were found on the S3 chip in both the 2010 obser¬ 
vation and the corresponding new observation. These sources 
are highlighted in white and the pulsar is labeled as P in Fig¬ 
ure [j] 

Careful consideration and modeling of the point-spread 
function (PSF) must be undertaken to reduce the effect of 
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Figure 2. Offset in the position of the pulsar in RA {top) and DEC {bottom) 
in the new observations from the position of the pulsar in the archival obser¬ 
vation plotted against the time since the first observation of PSRJ1741 with 
Chandra (years). The dashed line corresponds to the line of best fit in which 
the slope corresponds to the proper motion of the pulsar. 


Table 1 

Chandra observations of PSR J1741-2054 


Observation ID Exposure Time (ks) Observation Date 


archival 



11251 

48.78 

2010-05-21 

new observations 



14695 

57.15 

2013-02-06 

14696 

54.30 

2013-02-19 

15542 

28.29 

2013-04-01 

15638 

29.36 

2013-04-02 

15543 

57.22 

2013-05-15 

15544 

55.73 

2013-07-09 


changes in the PSF shape on the count distribution of our 
point-sources. To improve the astrometric accuracy and re¬ 
duce this effect, we simulate a PSF for each point source po¬ 
sition, in each observation, and use it to fit for the position of 
the source. 

To simulate a PSF of a point source, we use the software 
suite SAOT ratvFl which is designed to simulate the propaga¬ 
tion of photons from astronomical objects through the optics 
of the Chandra X-ray satellite. We use the aspect solution file 
of each observation and provided spectral information for the 
ray trace by extracting the spectrum of each stellar source us¬ 
ing the CIAO task specextract. We model the spectra using an 

11 http://cxcoptics.cfa.harvard.edu/SAOTrace/Index.html 
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Figure 3. Deconvolved image of J1741-2054 from ObsID 11251 obtained 
using arestore after 50 iterations. The cyan region corresponds to the PSF 
asymmetry seen in Chandra data when pushing to sub-ACIS-pixel resolution. 
The image has been logarithmically stretched and the colour enhanced to 
highlight the observed features. 


absorbed Mekal model in SHERPA. To improve the accuracy 
of the PSF modeling, we increase the normalisation of this 
spectrum by a factor of 100 before passing it into SAOT race. 
A model of a point source at its position is obtained by pass¬ 
ing the raytrace from SAOT race into the program ;V/,4 //aQ 
Each PSF model is corrected for the science instrument mod¬ 
ule (SIM) offset from nominal location and filtered using the 
Good Time Interval (GTI) data from the original event file. 

To determine the position of the sources we use the maxi¬ 
mum likelihood “figure of merit” (FoM) technique developed 
bv iVan Etten et all (120121) . We generate a 39 pixel by 39 pixel 
image of the modeled PSF, binned to 1/9 ACIS pixel resolu¬ 
tion. This is then compared to 0.3 - 5.0 keV source images of 
the same size but binned to native ACIS pixel resolution. PSF 
models and source images are produced for all observations 
and for each registration source. The PSF is shifted along the 
x and y axes of the 1/9 pixel grid and rebinned to native ACIS 
pixels. We then compute the FoM at each offset in pixel co¬ 
ordinates, giving us a map of the likelihood of the observed 
counts with respect to the x and y position. To determine the 
best-fit position of the source, we fit a two dimensional Gaus¬ 
sian to the FoM surface, with the minimum of this surface 
providing the best fit position of the source. The standard de¬ 
viations along x and y are estimated by calculating the square 
root of the eigenvalues of the covariance matrix. 

Prior to registration we checked each source for any opti¬ 
cal counterparts using ViziePFl. Sources 6, 8, 9, 10 and 13 
all have optical counterparts within 2" of the PSF-fit position. 
Source 8 and 13 have estimated proper motions of (pra,Pdec) 
= (22,-62) mas yr" 1 and (20,-22) mas yr" 1 , respectively (er¬ 
rors Sp.RA ~ 7 mas yr -1 , SpoEC ~15 mas yr" 1 ). We corrected 

12 http://space.mit.edu/CXC/MARX/ 

13 http://vizier.u-strasbg.fr/viz-bin/VizieR 
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Figure 4. Logarithmically scaled, merged Chandra exposure-corrected 0.3-5.0 keV ACIS-S image of the pulsar and its extended emission. Left : The green 
circular region defines the 0.70 // region used to extract the spectrum of the pulsar from all available observations. The cyan rectangular region defines the 
region used to extract the spectrum of the pulsar’s compact nebula from all observations. Right : The green regions define the apertures used to look for spectral 
variability in the extended emission of the pulsar. The pulsar was excluded from the spectral extraction of region 1. 


Table 2 

Frame shifts and their uncertainties used for registration. 


ObsID 

RA 

(arcsec) 

DEC 

(arcsec) 

14695 

-0.02 ± 0.02 

0.48 ± 0.02 

14696 

-0.01 ± 0.02 

0.25 ± 0.02 

15542 

0.02 ± 0.04 

0.43 ± 0.04 

15638 

-0.06 ± 0.03 

0.36 ± 0.03 

15543 

0.23 ± 0.02 

0.29 ± 0.03 

15544 

-0.11 ±0.03 

-0.29 ± 0.03 


for these nominal proper motions, but also confirmed that our 
final astrometric solution was insignificantly changed if we 
excluded these two stars from the analysis. 

Using the FoM positions of our registration sources and 
their uncertainties as a reference grid to perform the relative 
astrometry, we determined the best translation transformation 
needed to register the images from the two epochs using the 
CIAO command wcs_match. For each observation we adopt 
the Chanel ra-determined roll angle. In Table [2] we list the 
best-fit frame shifts and their uncertainties. The uncertainties 
in these shifts were calculated by adding in quadrature the er¬ 
rors in the differences between each source before and after 
shift. Adding a rotation to the transformation did not produce 
a statistically significant improvement on the best-fit transla¬ 
tion. 

To calculate the position of the pulsar after registration, we 
first calculate the position of the pulsar in the unregistered 
frames. We simulate a PSF at the position of the pulsar in 
each observation using the method described for the registra¬ 
tion sources. To define the energy dependence of the PSF, we 
extract a spectrum of the pulsar and fit it with an absorbed 
power law plus blackbody model (see Section 5). For each 
observation, we generate a 6 pixel by 6 pixel image of the 
PSF model of the pulsar that is binned to 1/4 ACIS pixel res¬ 
olution. Using an image of the 0.3 - 5.0 keV pulsar events of 
the same size and binning, we fit for the position of the pulsar 
in SHERPA. In this fit, we used a delta function (the pulsar) 
plus a two dimensional Gaussian (the circumpulsar PWN), 
both convolved with the PSF. The resulting pulsar fit position 


for each frame was then registered by applying the best-fit 
transformations (Table |2). 

To quantify the proper motion of the pulsar, we plot in 
Figure [2] the offset in the position of the pulsar between the 
archival observation and the new observations against the 
number of years since the archival observation. The uncer¬ 
tainties in the offset are calculated by adding in quadrature the 
uncertainties in the fit positions of the pulsar, the uncertainty 
in the frame shifts and the systematic uncertainty associated 
with choosing a particular tangent plane when creating an im¬ 
age in sky coordinates. We fit the offset using a linear function 
that corresponds to the positional shift of the pulsar between 
the archival observation and the new observations and this is 
seen as the dashed line in Figure 12 

We obtain a proper motion of [1 racos( 6) = -63 ± 12 
mas yr~ and /idec = 89 ±9 mas yr 1 . This corresponds to 
a total proper motion of 109±10 mas yr -1 . Assuming a dis¬ 
tance of 0.38 kpc to the pulsar, this translates to a transverse 
velocity of (196 ± 18)c/o,38 km s -1 . The position angle of the 
proper motion is 215° ± 6° east of north. The proper motion 
axis points in the opposite direction of the extended X-ray trail 
as expected (see Figure Q} 

4. X-RAY IMAGING 

Using the archival C handr a observation of J1741 (ObsID: 
11251), iRomani et al.l (1201 Oh performed a PSF subtraction 
of the pulsar point source to look for any small-scale struc¬ 
ture surrounding the pulsar. They discovered that the region 
around the pulsar appears to be slightly extended and they 
associate this feature with the equatorial torus_of the PWN . 

Using the same data as analyzed here, I Karpova et al, ‘ ( 2014 ) 
searched for evidence of such structure by performing fits to a 
two-dimensional Gaussian convolved with PSF models gen¬ 
erated for each observation. They find no evidence for any 
small-scale extended features other than for a small emission 
feature associated with a known mirror artifac0 We have 
carried out a similar investigation using image deconvolution 

14 http://cxc.harvard.edu/ciao/caveats/psf_artifact.html 















Table 3 

Spectral fit values of the extended emission of the pulsar defined by the cyan coloured region in Figure® and the green regions defined in Figure[4j>. All 

uncertainties are one cr. 
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Region 

AfolxlO-Tcnr 2 * 

r 

Absorbed Flux 
10“ 14 erg cm -2 s -1 

Unabsorbed Flux 
10“ 14 erg cm -2 s -1 

X 2 / dof 

compact PWN 

1 

2 

3 

4 

5 

1 20^ m 

1 - ZU - 0.07 

1.60 ± 0.20 

1 07+O.I8 

50^ 

17 ^ 

1 63+ 012 
j7q+8 :1A 
_ 1 /u -0.11_ 

?56« 

54 + 8 : 1 ? 

asi" 

4 oc+ 0.28 

540+8:^ 

_ 34 U - 0.36 _ 

72 1 C+ 0.09 

- 0.07 

1.81T0.33 

1.673=0.16 

1.10±0.11 

4 01+ 0.04 

603 -O.I6 

0.90 

0.99 

0.84 

0.80 

0.99 

1.16 


* Fixed at value from joint compact nebula. 


techniques. We simulate a PSF of the pulsar in each obser¬ 
vation using MARX following re-reduction of our Chandra 
observations with CALDB 4.4.7 to match the calibration data 
used for MARX. We define the energy dependence of the PSF 
as described in Section 3 and use the dither pattern of the 
observation with an aspect blur of 0.07, which corresponds 
to the uncertainty in the telescope pointin^Q We also cor¬ 
rect for SIM offset. Using this PSF we deconvolve a 0.3 - 
5.0 keV pulsar image using the Lucy-Richardson deconvolu¬ 
tion algorithm dLucvl 19741) implemented using the Cl AO task 
arestore. All images were binned to quarter-ACIS pixel. We 
ran arestore multiple times using a number of different iter¬ 
ations between 10 and 200 to determine convergence of new 
features. No new structures appeared after 50 iterations. 

The deconvolved image from one observation is shown in 
Figure 0 The emission is well described by a point source, 
with the exception of the artifact feature (identified in cyan). 
We thus see no conclusive evidence of other small-scale struc¬ 
ture in the immediate region surrounding the pulsar, consis¬ 
tent with the results reported bv lKarpova et ah ' ( 2014 )). 

5. SPECTRAL ANALYSIS OF THE X-RAY EMISSION OF THE 
PULSAR AND ITS TRAIL 

5.1. Pulsar 

As noted above, modeling of the Chandra PSF for use in 
determining an accurate pulsar position requires knowledge 
of the source spectrum. XMM-Newton observations establish 
a two-component spectrum for the pul sar, w ith a blackbody 
accompanied by a power lawjMarelli etjil.1 12014 9. Similar 
results were derived bv lKarpova et al.l (120141) using the same 
Chandra data reported here. We have re-analyzed these data 
by extracting events from a 0.70" radius region centered on 
the pulsar, shown as the green circle in Figure |4ji. The size of 
this region was chosen to minimise the contamination from 
the PWN, and subsequent modeling results have been cor¬ 
rected for the finite encircled energy fraction. All spectra 
were grouped with a minimum of 20 counts per bin, and a 
background spectrum was obtained using a source-free circu¬ 
lar region with a radius of 30", as shown by the magenta circle 
in FigureQ] Using the CIAO pileupjnap, we determined that 
approximately 5% of the pulsar events suffer from pileup. We 
thus included a pileup model in our spectral fits, where the 
frame time is 1.7s and the PSF fraction is allowed to vary. All 
other parameters in the pileup model were frozen at default 
values. 

To constrain the column density, we fit spectra from the 
compact nebula to the northeast of the pulsar (see Section 5.2) 
using an absorbed power law (more details in section 5.2). 
We obtained Nh = (1. 20^07) x 10 21 cm -2 , in good agreement 

15 http:/ / space . mit. edu/CXC/marx/news . html 


with the above studies, and adopt this value in all our models 
of the pulsar spectrum. We use the Wilms et al. abundance 
table throughout our analysis (IWilms et al.ll2000 ). Modeling 
the spectrum with a power-law plus blackbody model, we ob¬ 
tained a photon index of F = 2.68 ± 0.04 and a blackbody 
spectrum with a temperature of kT et f = 61.7 ± 3.0 eV, in ex¬ 
cellent agr eement with th e resul ts reported by iMarelli et alJ 
(120141) and iKarpova et al.l (120141) . Omitting the pile-up cor¬ 
rection yields similar values. We use the best-fit values above 
for PSF modeling of the pulsar. 

We also ran fits using magnetised neutron star atmosphere 
models ( nsa and nsmax in SHERPA) for the thermal com¬ 
ponent. These gave somewhat different temperatures and 
emitting areas, but did not significantly improve the qual¬ 
ity of the fit. For example a magn etic carb on atmosphere 
model (nsmax model 12006, lMori & Ho (200 7)) gave a tem¬ 
perature of kT e fj = 86.0 ± 9.0 eV and emission radius of 
Remis = (4.90! 2 -°Mo.38 km. The power law component was 
only slightly affected with F = 2.63 ± 0.03. 

5.2. PWN and its extended X-ray trail 

To analyse the spectrum of the compact nebula described 
above, we extracted spectra from each observation using the 
cyan rectangular region in Figure and combined these us¬ 
ing specextract. We used the same background spectrum as 
for the pulsar spectrum. The compact nebula contains a total 
of ~900 counts and we binned the combined spectrum with 
a minimum of 20 counts per spectral bin. The PWN spec¬ 
trum is consistent with an absorbed power law with an index 
ofT= 1.60 ±0.20. 

To determine whether there is any spectral variation in the 
PWN and its extended trail, we extract spectra from the five 
regions defined in Figure Et*. These regions correspond to 
roughly the same spectral regions reported by IMarelli et alJ 
(120141) . except that we investigate smaller regions in the com¬ 
pact nebula near the pulsar. In the outer portions of the nebula, 
the count rate is too low to obtain good spectra in smaller re¬ 
gions. We model each region individually using an absorbed 
power law, where we fix column density to the value derived 
earlier from fits to the inner nebula region but let photon in¬ 
dex and the normalisation vary. We have listed in Table 0 the 
absorbed flux, unabsorbed flux and the reduced for eac h 
region, as well as the best-fit parameters from fitting the trail 
spectra. There is no evidence of systematic variation in the 
photon index of the compact nebula and trail (regions 1 - 5). 
The photon index of region 1 is slightly higher than that of 
the other four regions, but is consistent within uncertainties 
with all regions except for region 2. This slight variation be¬ 
tween region 1 and 2 could suggest that region 1 is affected 
by leakage of the softer emission from the point source. Mod¬ 
eling the photon index as a function of distance from the pul- 
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Figure 5. The past trajectory (red) of PSR J1741 —2051 in Galactic coordinates. Positions are marked every 10 5 y, with uncertainty ellipses at t c = 3.91 X 10 5 y 
and at 10 6 y. For comparison we show the locations and sizes of OB associations from the Mel'Nik & Efremov! ( 1995 ) catalog (blue dotted ellipses) and the 
Hipparcos sample of OB stars (green circles). A substantial concentration of OB stars at ~ 0.3 — 0.6kpc lies near the pulsar track at ~ r, and a more distant 
(d = 1.45 kpc) OB association overlaps with the track uncertainty at this age. 


sar using a linear regression fit with a constant function in 
SHERPA, we obtain T = 1.67 ±0.06 for the trail. The global 
PWN index derived from the SHERPA fit and the valu es w e 
obt ained in Table |3l are c onsistent with iKarpova et all (120141) 
and iMarelli et al.l 720141) . who derived T = 1.74 ±0.07 and 
T = 1.78 ±0.15 respectively for the PWN. 

6. DISCUSSION 

Using Chandra observations of PSR J1741-2054 that span 
~ 3.2 year period, relative astrometry measurements have 
identified a proper motion of/t=109±10mas yr“ 1 . This corre¬ 
sponds to a modest velocity of (196 ± 1 8)t/n k m s" 1 , which 
agrees well with the velocity derived by iRomam^ahy gOlOj) 
using Her spectroscopy. The larger distance of IKarpova et al.l 
(2014 ) gives a transverse velocity of ~ 400 km/s, inconsistent 
with that obtained from optical spectroscopy. The direction 
of the proper motion i s 205°±6° east of nor t h, opposite the 
elonga ted X-ray trail. iNg & Roman i (2004), iJohnston et al.l 
(120051) and lNg & Romanil (120071) found that the direction of 
proper motion of a pulsar is, generally, approximately paral¬ 
lel to its rotation axis. 

In Figure [5] we plot the pulsar track (red) in Galactic co¬ 
ordinates (note the expanded b scale). This has the pulsar 
skimming above the plane. It does not intersect the plane it¬ 
self unless one extrapolates an unreasonable ~2x 10 6 y; how¬ 
ever, the track starts within the ~ 50 pc OB star scale height 
for t ~ t c and distances d < 1 kpc. For comparison we plot 
the positions of the Hipparcos catalog OB stars (green), with 
circle size proportional to the parallax. This set is quite com¬ 
plete, with useful parallaxes, to ~ 500pc, and increasingly in¬ 
complete at large r distances. At large distan ce the cataloged 
OB associations (iMel’Nik & Efremov l l l 995 ) provide plausi¬ 
ble pulsar birthsites, and their cataloged extent is plotted by 
the blue dotted ellipses. Intriguingly, one association over¬ 
laps the pulsar track, but this is at a likely unreasonable cat¬ 
alog distance of 1.45 kpc. We conclude that, with the pulsar 
motion passing along the Galactic plane, there will be many 
superposed massive star locations, and no definitive birthsite 


can be identified. However, there are certainly many plau¬ 
sibly associated massive stars consistent with our preferred 
d ~ 0.4 kpc, especially considering that some pulsar progeni¬ 
tors may be OB runaways with significant offset during their 
pre-explosion lifetime. 

Neither our deconvolved nor our PSF-subtracted images in¬ 
dicate conclusive evidence of small-scale structure surround¬ 
ing the pulsar that might be associated with a torus o r jet¬ 
like feature. The equatorial torus structure that lRomani et al.l 
( 2010 ) associate with a diagonal excess seen ~ 0.75" from 
the core of the puls ar image seems to have arisen from the 
mirror asymmetry. IKarpova et al. (120141) perform a similar 
analysis and come to the same conclusion. However, it is in¬ 
teresting to compare the nebula head and proper motion with 
the Ha structure described in Romani et al (2010). In Fig¬ 
ure [6] we see that the pulsar lies very close to the bow shock 
limb (accuracy limited by our relative X-ray/optical astrome¬ 
try). Interestingly, our measured proper motion is consistent 
with, although nominally slightly south of, the Ha nebula’s 
symmetry axis. However the X-ray PWN trail fills only the 
southern half of the apparent Ha- cavity, punching out through 
a gap at the back end of the Ha emission and continuing to 
the arcmin-scale trail beyond. The origin of this asymmetry is 
unclear, but a clue may be seen in the X-ray contours, whose 
ridge line lies at PA« 70°, i.e., misaligned with the proper mo¬ 
tion by (=s 35°. This suggests a second symmetry axis in the 
PWN, possibly due to a pulsar jet or other outflow concentra¬ 
tion. This directs the shocked PWN plasma to the southeast, 
preferentially filling this half of the Ha cavity. A more com¬ 
plete discussion of the PWN geometry, including the 3D Ha 
kinematics, is in preparation. 

The X-ray spectrum of the pulsar requires a combination 
of non-thermal and thermal model components. The emis¬ 
sion is dominated by the non-thermal component (~ 75% of 
unabsorbed flux), indicating that the majority of the X-ray 
emission is magnetospheric in nature. The emitting radius 
implied by the blackbody model for 11741 corresponds to 
( 4 . 5^2 5)^038 km. This is substantially smaller than any vi- 
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Figure 6. The head of the PSR J 1741—2051 nebula. The red channel shows 
Hen emission at 6562-6564A while the green channel shows average of the 
continuum 6551 —6561,6565—6575A with data abstracted from an AAT SPI¬ 
RAL IFU observation. The blue channel show the 0.3-7keV X-ray photons 
from our combined new observation. Green contours help one to read the 
X-ray structure, while the yellow arrows show our measured proper motion, 
with the arrow heads marking the position in 100 yr. 

able neutron star radius dLattimer & Prakashi : 2007 ). suggest¬ 
ing that this thermal emission arises from ho t spots o n the sur¬ 
face, plausibly near th e ma gnetic poles dHo & Heinkell2009l) . 
In fact, iMarelli et al.l (120141) do detect a pulsed thermal com¬ 
ponent for PSR J1741, also supporting such a surface temper¬ 
ature inhomogeneity. 

The X-ray emission from the compact nebula and the trail 
is consistent with an absorbed power law. There is no dis¬ 
cernible evidence of spectral variation with distance from the 
pulsar and the spectrum of the entire tail can be described 
by T = 1.67 ±0.06. We compute the minimum (equiparti- 
tion) energy by approximating the X-ray emission from the 
PWN (region 2, 3, 4 and 5 in Figure 0J)) as a cylinder with 
length / ~108 arcsecs (0.20c/o.38 pc) and width w ~±8 arcsecs 
(0.03r/o.38 pc), comprising a volume of V ~5.0x 10 5 V^0.38 
cm 3 , where <j> is filling factor. The minimum energy in rel¬ 
ativistic particles and magnetic field required to p roduce a 
synchrotron source of a given luminosity dPacholczvkll 19701) 
yield E mm ~ C(1 + L A syn , where k is the ion to elec¬ 

tron energy ratio, L syn is the synchrotron luminosity and C 
is a function dependent on energy, electron charge, speed 
of light and the mass of the electron in Gaussian cgs units 
(see iPacholczvk] dl970i) f. In the following, we considering 
only the leptonic case, where n = 0. The total luminosity of 
the PWN is L(0.5-10.0 keV) = 2.36 x 10 30 erg s“\ giving 
E min ~ 5.50 x lO'V/ 7 ^ 7 erg. The associated minimum- 
energy magnetic field is ~ (D(l + K)L syn ) 2 / 1 V~ 2 ^ 1 , where 
D is a function similar to C. This magnetic field is ~ 
15(/r 2 / 7 r/ 0 38/ttG, leading to a lifetime of the X-ray emitting 
leptons of T syn ~6.4x 10 4 B^E k jy 2 years or ~ llOOyr at an 


observed photon energy of 1 keV. This is comparable to the 
(length/proper motion) = 108"/109 ~ 10 3 yr required for the 
pulsar to traverse the bright trail with our observed proper mo¬ 
tion. Thus it is not surprising that there is no dramatic spectral 
steeping along the trail. If the PWN electrons flow at even 
faster speeds within the trail, this conclusion is even stronger. 

7. CONCLUSION 

Using ~ 300 ks of Chandra ACIS-S observations of PSR 
J1741-2054, we were able to determine the proper motion of 
the pulsar with a detection significance > 3er. The direction 
of the proper motion is aligned with the extended PWN emis¬ 
sion, and corresponds well with a symmetry axis of the asso¬ 
ciated Her nebula. The diffuse X-ray emission immediately 
behind the pulsar is concentrated in the southeastern portion 
of the Ha nebula, possibly suggesting another flow axis from 
a jet or torus in the pulsar system. The trajectory of the pul¬ 
sar, extrapolated over the characteristic age, does not provide 
a compelling correlation with known OB associations at the 
distance of the pulsar, although there are many massive stars 
consistent with this distance that could potentially have had a 
common origin. 

The pulsar spectrum is well described by an absorbed 
power law accompanied by a blackbody with an emission ra¬ 
dius of (4.5 !^)Jo. 38 km and a temperature of kT e ff = 61.7 ± 
3.0 eV, as found in earlier works. The thermal component, 
a hot region on the neutron star surface, is augmented by a 
magnetospheric or unresolved PWN power law component. 
The PWN plus its extended trail can be well described us¬ 
ing an absorbed power law and there is no evidence of vari¬ 
ation in the photon index with distance from the pulsar. The 
integrated luminosity of the PWN over the 0.5 - 10 keV is 
2.36 x 10 30 erg s -1 . This represents 0.02% of the pulsar spin 
down power, which is not atypical. We find no conclusive ev¬ 
idence of small-scale structure surrounding the pulsar that we 
can associate with a torus or jet-like structure. 
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